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The spatiotemporal response of crystals in x-ray Bragg diffraction resulting from excitation by an 
ultra-short, laterally confined x-ray pulse is studied theoretically. The theory presents an extension of 
the analysis in symmetric reflection geometry [JJ to the generic case, which includes Bragg diffraction 
both in reflection (Bragg) and transmission (Laue) asymmetric scattering geometries. 

The spatiotemporal response is presented as a product of a crystal-intrinsic plane wave spatiotem- 
poral response function and an envelope function defined by the crystal-independent transverse pro- 
file of the incident beam and the scattering geometry. The diffracted wavefields exhibit amplitude 
modulation perpendicular to the propagation direction due to both angular dispersion and the dis- 
persion due to Bragg's law. The characteristic measure of the spatiotemporal response is expressed 
in terms of a few parameters: the extinction length, crystal thickness, Bragg angle, asymmetry 
angle, and the speed of light. 

Applications to self-seeding of hard x-ray free electron lasers are discussed, with particular em- 
phasis on the relative advantages of using either the Bragg or Laue scattering geometries. Intensity 
front inclination in asymmetric diffraction can be used to make snapshots of ultra-fast processes 
with femtosecond resolution. 

PACS numbers: 41.50.+h,42.25.-p, 61.05.cp, 07.85.Nc 



I. INTRODUCTION 

The spatiotemporal response from crystals in symmet- 
ric x-ray Bragg diffraction in reflection (Bragg) geome- 
try has been studied in our recent publication [1 . Here, 
we extend the analysis to the generic case of asymmetric 
Bragg diffraction, both in reflection (Bragg) and trans- 
mission (Laue) geometries. "Asymmetric" means that 
the Bragg reflecting atomic planes are not parallel to the 
crystal surface. 

Understanding the time dependence of x-ray Bragg 
diffraction in crystals has attracted much attention since 
the 1990s, as the advent of ultra-fast (femtosecond short) 
x-rays pulses become a close reality. The temporal and 
spatial dependence of diffraction was first calculated us- 
ing the time dependent Takagi-Taupin equations. In par- 
ticular, an analytical solution for the Bragg reflected 
wave from an infinitely thick crystal was derived by 
Chukhovskii and Forster [3J. Numeric calculations of 
the time dependence of Bragg diffraction from a crys- 
tal heated by a laser pulse was performed in [3]. Cal- 
culations of the time dependence by Fourier transform- 
ing the known monochromatic plane- wave solutions from 
the classical dynamical theory |iHTTj have been consid- 
ered in several publications [T2lR9| . In particular, Shas- 
tri et al. [12, 13J performed numerical calculations, which 
have revealed signature features of time dependences of 
Bragg diffraction from crystals both in the Bragg-case 
and in the Laue-case geometries. Graeff and Malgrange 
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p~5l H6] obtained analytical solutions for the time depen- 
dence of Bragg diffraction in the Laue geometry, with the 
refraction effects at the crystal exit surface taken into ac- 
count. Bushuev [30] used Fourier transformation of the 
plane-wave solutions both in the frequency and momen- 
tum space, with the second order corrections included 
to more accurately account for the refraction effects, to 
obtain solutions in time and space and analyzed specific 
cases using numeric calculations. 

The present paper is focused on the development of the 
theory and on the analysis of the spatiotemporal response 
of crystals in Bragg diffraction to the excitation by an 
ultra short in time and spatially confined x-ray pulse in 
the general case of asymmetric reflection (Bragg) and 
transmission (Laue) scattering geometries. The primary 
goal of the present study is to understand the general 
phenomenology of the spatiotemporal response by un- 
covering the dominant underlying physics and identify- 
ing the key physical parameters that determine the char- 
acteristic time and space scales involved. For this pur- 
pose, we derive comprehensive solutions that can be writ- 
ten in the general case as a product of two independent 
envelope functions: the first is a spatiotemporal plane- 
wave response function that depends only on the crystal 
and scattering geometry, while the second is an enve- 
lope that is specific to the initial conditions of the inci- 
dent field. We derive analytical solutions for the response 
functions under several representative conditions, which 
clearly identifies the key physical parameters and makes 
possible a relatively simple interpretation of the general 
solution. 

When an ultra-short x-ray pulse instantaneously excites 
a perfect crystal, the output field is delayed and spread in 
time. The underlying reason behind this phenomenon is 
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that each frequency component excites a monochromatic 
eigen-wavefield in the crystal that propagates along its 
own direction with its associated group velocity. As a 
consequence, the time response is intrinsically connected 
to the lateral spatial distribution of x-rays leaving the 
crystal upon Bragg diffraction both in the reflection or 
forward directions. 

The paper is organized as follows. Comprehensive so- 
lutions for the spatiotemporal dependences of Bragg 
diffraction in reflection (Bragg) and transmission (Laue) 
asymmetric geometries are derived in Sec. [n] In par- 
ticular, in Sec. |II B| the solutions are derived for inci- 
dent ultra-short x-ray pulses with an unbounded plane 
wave front, and in Sec. |II C| for incident ultra-short x-ray 
pulses with a bounded wavefront. The solution for the 
bounded wavefront is a product of a crystal-intrinsic and 
geometry-specific plane wave spatiotemporal response 
function and an envelope function defined by the crystal- 
independent transverse profile of the incident beam and 
the scattering geometry. The response functions in the 
asymmetric Bragg geometry are derived analytically in 



Sect. Ill in the approximation of a non-absorbing and 
thick d ^> A'h crystal (the crystal thickness d being much 
larger than the characteristic extinction length of Bragg 
diffraction A# to be more precisely defined below). The 
response functions in the asymmetric Laue geometry are 
derived in analytical form in Sect. |IV| for a non-absorbing 
crystal of arbitrary thickness. Applications of the theory 
for self-seeding of XFELs and for ultra-fast time mea- 
surements are discussed in Sec. IVl 



II. COMPREHENSIVE SOLUTIONS FOR 
SPATIOTEMPORAL CRYSTAL RESPONSE 



We study here the spatiotemporal dependence of Bragg 
diffraction of ultra-short, laterally bound x-ray pulses 
from a system of parallel atomic planes in a flat crystal 
plate. Generic solutions are derived in three consecutive 
steps. First, well-known solutions of the dynamical the- 
ory of x-ray Bragg diffraction in crystals [3HTTI 1TB] for 
incident monochromatic plane waves are briefly summa- 
rized in Sec. ( II A ) . In Sec. II B we derive solutions for an 



initially ultra-short incident pulse with boundless plane 
wavefront. Finally, solutions are obtained in Sec. |II C| for 
an ultra-short incident pulse with confined wavefront. 



Monochromatic Plane Wave Solutions 



One of most fundamental results of the dynamical theory 
of x-ray diffraction in perfect crystals is the concept of 
monochromatic eigen wavefields in crystals introduced by 
Ewald almost 100 years ago A similar concept in the 
electron theory of solids was introduced later by Bloch 



in 1928, which are therefore generally known as Bloch 
waves. In the simplest case, an incident monochromatic 
plane wave £ i exp(iK n i — iw i), with a frequency ui and 
wavevector K = (lo /c)u , propagating along the opti- 
cal axis tt , excites in the crystal a waveficld 

2>(r, t) = £ 1 exp(-kV) V R 0H exp [i(K + H)r] (1) 



which is a sum of plane wave components with wave vec- 
tors K a + H. In the general case this sum involves all 
the reciprocal crystal lattice vectors H of the crystal. 
For each H there is a set of parallel atomic planes in 
the crystal perpendicular to H with an interplanar dis- 
tance d H = 2n/H which actually composes the grating 
on which x-rays diffract. 

In the following we will consider the so-called two-wave 
case, where only two plane wavefield components are 
taken into account: the wave associated with forward 
Bragg diffraction H = 0, and one Bragg diffraction com- 
ponent with nonzero H, for which \K a + H\ ~ \K \ = 
K a , and for which therefore the relative difference 
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has a very small magnitude. In particular, if a = we 
obtain Bragg's law 2K H + H 2 = 0, which can be also 
written as 



2,?Csm6> = H. 



(3) 



Here, 9 is the glancing angle of incidence to the atomic 
planes, which equals the angle between K and the 
atomic planes such that K H — — K Q H sin#. The quan- 
tity a ([2]) is an important parameter of the theory known 
as the deviation parameter, since it represents the devi- 
ation from Bragg's law. 

We take the atomic planes associated with the recipro- 
cal crystal lattice vector H to be oriented at an arbi- 
trary (asymmetry) angle 77 with respect to the crystal 
surface, as shown schematically by white parallel lines in 
Figs. [T|a) and [TJb) . Figure [TJa) represents diffraction 
in the reflection or Bragg scattering geometry, which is 
characterized by the diffracted wavefield exiting the crys- 
tal on the same side as the incident wave, while the for- 
ward diffracted wavefield propagates along the incident 
wave direction and exits the crystal on the opposite side. 
Figure [ljb) represents diffraction in the transmission or 
Laue scattering geometry, for which both the diffracted 
and forward diffracted wavefields exit the crystal from the 
surface opposite that of the incident wave. The crystal 
surfaces are defined by the unit normal vector z internal 
to the entrance surface. We restrict our theory to the 
case in which K a is directed in the plane composed by 
H and z, hereafter referred to as dispersion plane. 

The dynamical theory of x-ray diffraction allows one to 
calculate for each incident monochromatic plane wave 



3 




FIG. 1: Schematic presentation of two-beam x-ray Bragg diffraction from crystals (a) in the reflection (or Bragg) scattering 
geometry, and (b) in the transmission (or Laue) scattering geometry. The glancing angle of incidence to the reflecting atomic 
planes is 9, and the angle between the reflecting planes and the crystal surface (the asymmetry angle) is rj. The propagation 
direction u H of the Bragg reflected beam composes an angle i\> H with the internal normal z to the crystal surface. The angle 
■00 between z and the direction u of propagation of the incident beam is defined by the relationship %p + ip H = 29. For the 
scattering configurations shown in (a) and (b) ip H = n /2 + 9 — r], and therefore ip = 9 + -q — n/2. The permitted range of the 
asymmetry angle r\ in the Bragg-case geometry is — 9 < rj < 9, while in the Laue-case geometry it is 9 < r\ < tt — 9. Negative rj 
in the Bragg-case geometry corresponds to a configuration complimentary to that shown in (a) with the incident and reflected 
beams reversed. Pink and light blue areas indicate regions where the forward diffracted and diffracted beams can propagate. 
See text for other details and definitions. 



component exp{— i[ui t — K r]} both the monochromatic 
wavefield of forward Bragg diffraction 

2>( m V ,i) = E x ~ K r a ] RM (4) 

at any point r on the rear surface of the crystal, and the 
monochromatic wavefield of Bragg diffraction 

2>W(r H> t) = fi.e-H*- (*o +H)r H ] RqM _ (5) 

The Bragg diffracted field §5§ is determined at any point 
on the entrance surface r H in the case of Bragg geometry 
shown in Fig.Qa), or at any point on the rear surface r H 
for the Laue geometry shown in Fig.JlJb). Here, R on (u) ) 
are diffraction {T-L — H) or forward diffraction {Ti = 0) 
crystal amplitudes which are functions not only of ui a 
(whose dependence we explicitly indicate since it is most 
relevant to the discussion of the time behavior), but are 
also functions of the crystal thickness d, the direction and 
magnitude of K , the asymmetry angle rj, etc. At this 
point we are concerned with deriving general expressions 
for the spatiotemporal response of the crystal in x-ray 
Bragg diffraction, and so do not yet specify these reflec- 
tion amplitudes; explicit expressions for the amplitudes 



The in-crystal monochromatic wave field components 
T>^ (r n , t) given by Eqs. Q and ^ can also be used to 
calculate the field at any point r outside of the crystal. 
Using the continuity of the wave fields at the crystal- 
vacuum interface determined by the extremities of the 
vectors r u , the forward diffracted and diffracted wave 
fields in an arbitrary point r in vacuum can be written 
as 



(r,t) = T>W(r H ,t) 



(6) 



where K n is the wavevector of the forward diffracted 
{T-L = 0) or diffracted {Ti = H) field in vacuum. To 
match phase fronts the components of the in-crystal wave 
vectors can differ from the vacuum wave vectors only by a 
component along the crystal normal z. Since we assume 
that the crystal entrance and exit surfaces are parallel, 
this component is zero for the wavevector K Q , and makes 
it equivalent to the vacuum wavevector K a of the incident 
plane wave. However, this component is not zero for the 
vacuum wave vector of the diffracted wave. In the general 
case it can be written as 



K H =K +H + A H z. 



(7) 



R 0Ji are presented in Sec. Ill and Sec. IV 



The component A H z can be understood as an additional 



4 



momentum transfer due to refraction at the crystal vac- 
uum interface, and Eq. ^ as momentum conservation 
in scattering from the crystal. Since Bragg diffraction 
is an elastic scattering process and the vacuum is homo- 
geneous, the magnitude of the vacuum wavevector K H 
of the diffracted wave should be equal to the the vac- 
uum wavevector of the incident plane wave: \K \ — 
\K a \ = K a . From this condition and Eqs. j7||, Q we 
find: A H = K Q (—^ H ± yf*(* — a) (see [H] for details), 
where ^ H = z{K a + H)/K . For small a, the additional 



momentum transfer can be closely approximated by a 
Taylor expansion in a: 



A, 



K, 



(8) 



Now from Eqs. Q, j5|, ([6]), and Q we can write for the 
monochromatic forward diffracted and diffracted wave- 
fields 



S^\r,t) = e.e-^ot ~ (K a + H)r] e iA H (r - r H )z r^). 



(9) 
(10) 



Here (r — r H )z is the shortest distance from the obser- 
vation point r to the crystal surface. Since z is perpen- 
dicular to the crystal surface, it is actually independent 
of r H1 and (r — r H )z = rz if the extremity of r H is on 
the entrance surface, and (r — r H )z = rz — d if the ex- 
tremity of r H is on the rear surface. Here d is the crystal 
thickness. 



B. Ultra-Short Incident Pulse with Boundless 
Plane Wavefront 



To study the spatiotemporal dependence of x-ray diffrac- 
tion, in the next step we investigate the response to an 
initially ultra-short (instantaneous) x-ray pulse. We as- 
sume that the x-ray pulse propagates along the direction 
of the unit vector u which is in the dispersion plane built 
by vectors H and z, and that the propagation direction 
u a makes a glancing angle of incidence 9 with respect to 
the reflecting atomic planes. 

The x-ray pulse is ultra-short, has a vector amplitude 
5. , and has an infinite extent in the transverse direction 
i, 1 ti,. In this case, the pulse at time t and spatial 
point r can be presented by the delta function <5(r) of 
the argument r = t — u r /c, with c the speed of light 
in vacuum. The latter is equivalent to an infinite sum of 
monochromatic plane wave components given by 



S(t) 



an 

2tt 



-i(cj + 0)r 



(11) 



T = t- 



u n r 



n = 



(12) 



Here, we single out one plane wave component with a fre- 
quency u (we assume to ^> CI), which we define to satisfy 



the condition a = defined in Eq. |2]). In other words, 
we are selecting out the frequency for which Bragg's law 
Q is fulfilled. For the frequency ui Bragg's law reads 
ujsaxO = Hc/2. With this convention, the deviation pa- 
rameter a ([2| can be presented as 



a = —4 — sin 



1 - 2- 
u 



and the additional momentum transfer A, 
A H _ 2 sin 2 6 fl 

LO/C 7 H LO 





sin 2 9 \ 









(13) 



. (14) 



where 



7^ 

1h 



7o 



7h 



H 

uj/c 



(15) 



are the so called asymmetry factor b, and direction 
cosines j = cos tp , ~f H = cos ip H . In almost all of what 
follows we retain only the terms linear in the small quan- 



tity Cl/ui for the expressions for a (13) and A H (14), but 
we also present a brief description of the physics of the 
quadratic terms and how they can be included. 



Linear approximation 



Time t — is defined hereafter as the moment when the 
wavefront hits the point r = on the crystal. Similar to 
( [TT| ) , the spatiotemporal response of the crystal in Bragg 
diffraction £ n (r,t), both for diffracted (T~L — H) and for- 
ward diffracted (H = 0) components, can be calculated 
as an integral (in fact, a Fourier integral) 



t) 



2tt 



£i m) (r,<) 



(16) 
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over the monochromatic components (|9|-([l0|). A similar 
procedure was also applied in the previous publications 
[THH]. Using ^- jio), we obtain 



£ H (r,t) = S iC -^ T nG 0H (U, 



where 



OO 1Q 



2tt 



c 

_^_2sin 2 ^ 



H 

c 



C7a 



(17) 
(18) 

(19) 
(20) 



Here, <5 KJf is the Kronecker delta, which equals one only 
if H = H , otherwise it is zero. 

The plane wave crystal response functions G 0H (£ w ) in 
( 18 ) represent the spatiotemporal dependence of Bragg 



diffraction (H — H) or forward Bragg diffraction (H = 0) 
to the excitation by a 5-function-short incident radiation 
pulse with boundless transverse wavcfront. 



Two spatiotemporal variables are introduced in (17) 



p0| ). The variable £ n (H = 0, H) Q is the argument of 
the response function (18). The spatiotemporal variable 
T n (H = 0, H) ( 19 ) is in the argument of the exponential 
function of Eq (17). 



The spatiotemporal variable t u (H — 0, H) represents 
the difference between the absolute time t and the time 
u H r /c the plane wavefront, propagating along the optical 
axis u H from r — 0, would need to reach an arbitrary 
point r outside of the crystal, assuming the propagation 
is in vacuum and the wavefront is perpendicular to the 
optical axis u H . In other words, t u is the time delay 
for the radiation field at point r we are interested in, 
as compared to the trivial propagation of the pulse in 
vacuum along u H . 

The variable £ is equivalent to r , while the meaning 
of £ H may not be immediately evident. To gain more 
insight, we rewrite £ H in an equivalent form using (19)- 
and Bragg's law u = flc/(2sin0) 



T, 



t h +2 sin 9 
2 d sin 2 9 



H 



sin 9 

1h ' 



Z5 F 



c\1h\ 

in Bragg geometry 

1 in Laue geometry. 



(21) 



The parameter T d is a characteristic measure of time in 
Bragg diffraction associated with the crystal thickness d. 
In the Laue-case geometry, T d is equal to the total dura- 
tion of forward Bragg diffraction, which is given by the 
difference in path lengths for the wave to propagate along 
0B A' and OA, as shown schematically in Fig. [2j The to- 
tal duration of Bragg diffraction in Laue-case geometry 



FIG. 2: Schematic presentation of plane wavefront paths 
(solid vector lines) in Laue-case Bragg diffraction and forward 
diffraction from a crystal plate. The total duration of forward 
diffraction (0B A' - 0A)/c = T d pi), while total duration of 
diffraction is (0AB' — 0B)/c — T d /b. The lateral spread of for- 
ward Bragg diffraction is A A' = dcos2#/7 ff , and the lateral 
spread of Bragg diffraction is BB' = dcos20/y . 




FIG. 3: Schematic presentation of plane wavefront paths 
(solid vector lines) with possible multiple internal reflections 
from the rear and front crystal surfaces in Bragg-case Bragg 
diffraction and forward Bragg diffraction from a crystal plate. 
The forward diffracted wave 0ABC originating from the front 
surface reflection in B is delayed by (0AB-0B")/c = T d |2l] , 
as compared to the primary forward diffracted wave 0B" . The 
diffracted wave 0AB originating from the rear surface reflec- 
tion in A is delayed by (0AB — 0B')/c = T d /\b\ compared to 
the primary diffracted wave OS'. 



is determined by the difference in path lengths 0AB' and 
0B, which equals T d /b as can be derived from schematic 
in Fig. [2] Although the last term in (21) vanishes in 
the Bragg-case geometry, the parameter T d continues to 
play an important role. Unlike the Laue-case geometry, 
diffraction in Bragg-case geometry is not limited in time, 
because multiple reflections from the front and rear sur- 
faces take place [3 [TUl EJ HI] as shown schematically in 
Fig. [3j Accordingly, the parameter T d is a characteris- 
tic measure of time associated with crystal thickness in 
the Bragg-case geometry, where it measures the time be- 
tween multiple reflections as explained in the caption to 
Fig.0 
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Using Fig. [I] we can express z and H/H in (21) in 
terms of the unit vector u H along the diffraction optical 
axis, and the unit vector v H perpendicular to the axis 
as follows: z — cosip H u H + smip H v H , and H/H = 
smOu H — cosdv H . With these, the third term in (21) 
can be presented as 



2sin6» 

D = 



H 



2 sin 9 sin r\ 
sm(9 — rj) 



= Dv H , 

= - (1 + &)tan( 



(22) 



Eq. (21) for t; H thus can be now written as 



v„r 



(23) 



where t h = t — u H r/c ( 19 1 



The quantity D in (22)- (23 1 is the normalized angu- 
lar dispersion rate. It is a measure of the variation 
of the propagation direction K H /K of the diffracted 
wave (Q as a function of the incident photon energy 
hK /c, assuming a fixed direction K / K of the incident 
wave vector. Indeed, using (|7|) we obtain 5{K H /K ) — 
—Dv H (8K /K g ), see [TS] for details. We note that the 
normalized angular dispersion rate D is zero only in 
Bragg-case symmetric scattering geometry with rj = 0. 
In all other cases, including the "symmetric" Laue ge- 
ometry (77 = 90°), it is nonzero. Depending on the sign 
of 1], D can take positive or negative values in the Bragg- 
case geometry. In the Laue-case geometry 7/ > and D 
is therefore always negative. 



From Eq. ( 23 1 we conclude that the spatiotemporal vari- 
able £ H = t h +Dv H r /c+T d contains in addition to t h an 
important term Dv H r/c which describes a spatially lat- 
eral (perpendicular to the diffraction axis u H ) amplitude 
modulation G 0H (£ H ) of the diffracted radiation field (18). 



The amplitude modulation occurs due to interference of 
different spectral components propagating in different di- 
rections, which arise from the angular dispersion due to 
the additional momentum transfer A H - Eq. Q and ( |T4| . 
The effect is very generic and vanishes only in one case, 
in symmetric Bragg geometry when 77 = - Fig [lja) . 
We will refer to this effect as the angular dispersive lat- 
eral spatial modulation of the diffracted wavefield (here 
by angle we mean the angle of reflection rather than the 
angle of incidence) . 



2. Nonlinear phase 



In our preceding discussion we have focused on the phase 
contributions that are linear in the frequency difference 
f2. While this is a very good approximation for the wave- 
fields in vacuum close to the crystal surface, the nonlinear 
(~ Q 2 ) contributions inherent in the additional momen- 
tum transfer A H give rise to additional physics over the 



potentially large propagation distances between the crys- 
tal surface and any experimental sample/ detect or. Here, 
we briefly quantify this effect and summarize its physical 
origin. 



Using Eq. |15|) , the definition of ip and ip H given in the 
caption of Fig. [I] along with the definition of D ( 22 ) , the 
expression for the additional momentum transfer ( 14 ) can 
be rewritten as 



A H _ 2 sin 2 9 O 
ui/c 7 H u 



D 2 



(24) 



Including the quadratic phase dependence in Eq. 
the wavefields (17)-(18) can be presented as 



(24), 



J-co 27T 



S(n) = exp 



LP 2 /? 
2cw 7ff ' 



a 2 



(25) 
(26) 

(27) 



where p = (r — r H )z is the shortest distance from the 
observation point to the crystal surface, so that p/j H is 
the propagation distance along the optical axis. We note 
that the expression for the nonlinear in CI phase factor 
S(Q) in (27) is in agreement with that obtained earlier 
by Bushuev in [20. , despite the different approaches used. 
With the help of the Fourier convolution theorem we ob- 
tain for X m (£J: 



poo 

Jo 



0. (28) 



F(tn " = 



W 2 p 



exp 



2p 



D 2 



(29) 



Here G on is the crystal response given by Eq. ( 18 1, while 



F(£) is the Fourier transform of S(Q). The convolution 
is similar in form to that associated with paraxial evolu- 
tion for the field G on (£), with the Fourier transform of S 
serving as the associated Green function, p/7 H the prop- 
agation distance along the optical axis, and £/D playing 
the role of the "transverse" coordinate. 

In fact, the factor 1/D gives the amount the inclined 
intensity front will spread in time due to natural vac- 
uum diffraction broadening along the transverse coordi- 
nate u H . In symmetric Bragg (reflection) geometry the 
reflected intensity and phase fronts are parallel, D = 
and the time structure remains invariant, while in all 
other cases D ^ and the separated, inclined intensity 
fronts will tend to smear together as the distance from the 
crystal p increases. The maximum distance over which 
the linear approximation holds and we can ignore this 
spreading can be estimated as 



D 2 oj 



U! 



(30) 
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where J7 B is the typical frequency range of Bragg diffrac- 
tion, i.e., the range in f2 over which R 0H (Q) is apprecia- 
ble. For example, if f2 B /w ~ 10 -4 , the radiation wave- 
length A = 2ttuj/c ~ lA, and j H /D 2 sa 1, the linear 
approximation breaks down at the rather small distance 
p R3 1 cm from the crystal. The smaller the bandwidth of 
the Bragg reflection f2 B , the larger is the distance from 
the crystal over which the linear approximation holds. 

In the following we will neglect these nonlinear effects due 
to vacuum diffraction upon propagation away from the 
crystal. If they have to be taken into account, one should 
replace the response function G 0H ( 18 ) with X 0H (£ u ) 
( 28 1- ( 29 1 in the equations presented below. 



The spatiotemporal response £ H (r,t) of the crystal in 
Bragg diffraction to the excitation by the ultra-short and 
laterally bound x-ray pulse ( 32 ) , both for diffracted (H = 
H) and forward diffracted (Ji = 0) components, can now 
be constructed as a Fourier integral over 9 of the plane 
wave solutions ( 17 )-( 18 ) : 



£u(r,t) =£ { 



t h (6) =t 



d9 
2^ 



u n (9)r 



t*0) =r H {6) + D{e)^ i ^+T d {0)5 E 



(33) 



Ultra-Short Incident Pulse with Confined 
Wavefront 



Here we are using again an important condition that the 
carrier frequency ui(0) satisfies Bragg's law 



In the next step, we introduce an incident x-ray pencil- 
beam directed along the unit vector u (8) = u . The 
wavefront of the pencil-beam is bounded in the direction 
v (9) = v perpendicular to u by the transverse profile 
II(i> ), where v — v r. We assume that the profile has a 
characteristic width of a v and can be written as a Fourier 
transform of the angular profile: 



n(« ) 



dQ 
2n 



T((9)cxp -w (u)/c){6-6) . (31) 



The characteristic angular spread a g in T(9) is related to 
cr v by the uncertainty relationship cr g a v ~ c/u>. In par- 
ticular, for a pencil-beam of x-rays with a photon energy 
hu ~ 12 keV (A = 2irc/ui ~ 0.1 nm) and a lateral spread 



of a v ~ 10 /xm, the angular spread a g 



10- 



rad. 



An ultra-short-in-time incident pencil-beam £ l (r,t) = 
£.e~ luJT o S(t )U(v ) can be presented as a Fourier inte- 



gral over 9 of the plane- wave components ( 11 1-( 12 1 prop- 



agating along directions u a {9) at glancing angles of in- 
cidence 9 to the atomic planes around the central angle 



d6 
2^ 



r(9)e-' luJT »^5 \t {9) 



T(0) 



T (fi)=t- 



3 ^ 

u (9)r 



2tt ' 



-i(« + n)r (fl) 



(32) 



The ultra-short-in-time pencil-beam presentation ( J32[ ) is 
valid provided the ^-dependence in the delta- function 
(5 t (9) can be neglected. Since t (9) ~ r — (6 — 9), 

this is valid if the time delays t we are considering are 
much longer than the inverse frequency: r 3> o g cr v /c ~ 
1/uj. 



Lj{9)sm9 = Hc/2, 
equivalent to ([3]), and a = condition. 



(34) 



Since only small 9 — 9 values are significant, we can use 
u u {9) = u u + 5u n , with 8u H = —v H 59. Here, the vec- 
tors v u are perpendicular to u n and directed as shown 
in Figs, [lja) and|ljb). Applying this result and Tay- 
lor expanding t u {9) and u)(6) from (33 1 to first order in 
\9 - 9\ < 1, we obtain 



r n (0) 



+ ^(9-8), 
c 

1 - (6- 0)cot( 



lo = oj{8), (35) 



t h (9)lu(8) —ojt h + — (v n — t h c cot ( 



G 0H (£ n ) is a slowly varying function compared to 
exp(— \ujt h ). Therefore, one can neglect dependence of 



G[£ M (9)] on 9 in performing integration over 8 in (33) 
provided the lateral shift v n or/and angular spread cr e 
are not too large, so that t w ^> a v /c. With these as- 
sumptions, and using again Eqs. (|2l[)-(23 1, we arrive at 



the following general expression for the spatiotemporal 
dependence of Bragg diffraction from a crystal, excited 
with an ultra-short-in-time pencil-beam with a lateral 
spatial distribution IL(v ): 



£ n (r,t) =£, G on {^ n )U(v H - t h c cot 0)e~ 



T =t — 



H = (0,H), 



(36) 



Equations (36) reveal an interesting general property: 



the spatiotemporal response in Bragg diffraction (H = 
H) or in Bragg forward diffraction {% = 0) is given by a 
product of the corresponding plane-wave spatiotemporal 
response function G 0H (£ H ) (|18[) and the spatiotemporal 
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envelope function H(v n — T n ccot8), whose peak shifts 
along v H perpendicular to the optical axis u H linearly in 
time. Thus, a fixed relationship (36) exists between the 



time delay of the crystal response and the peak of the 
lateral shift. In other words, the time delay is mapped 
onto the lateral spatial shift. 

In our previous paper [T], we have shown that the spa- 
tial shift takes place in symmetric Bragg diffraction in 



Bragg scattering geometry. The solution ( 36 ) generalizes 



that result to asymmetric diffraction both in reflection 
(Bragg) - Fig. [ija), and transmission (Laue) scattering 
geometries - Fig. [TJb). This result can be interpreted 
as follows: the incident wavefield with bounded wave- 



front is presented in Eq. ( 32 ) as a superposition of plane 
waves propagating at different angles of incidence. At a 
different angle, Bragg's law is fulfilled for different pho- 
ton frequency u{6) (34 1. As a result, the spatiotempo- 
ral response £ n (r,t) of the crystal in Bragg diffraction 
( 33 ) is a superposition of wave fields with different car- 



rier frequencies lo{9) (34 1, resulting in a lateral spatial 
modulation, or, equivalently, in a lateral spatial shift. 
An alternative interpretation of the derived above gen- 
eral relationship between the time delay and spatial shift 
is discussed in Appendix [Cj We will denote this effect 
as the lateral spatial modulation due to Bragg's law of 
dispersion, to distinguish it from the spatial modulation 
due to angular dispersion discussed in Sec. |IIB[ and refer 
to II() in Eq. (36) as Bragg's law dispersion envelope. 



We conclude: the spatiotemporal response of the crys- 
tal to the excitation with an ultra-short and laterally 
bounded x-ray pulse is accompanied by lateral spatial 
modulations driven by two different mechanisms: Bragg's 
law of dispersion and angular dispersion. We will illus- 
trate manifestation of these two mechanisms using par- 
ticular cases in Sec. Mil and Sec. IIVI 



Finally, using Eqs. (36) the spatiotemporal dependence 



of the intensity of Bragg diffraction from a crystal excited 
with an ultra-short-in-time pencil-beam having a lateral 
spatial distribution Tl(v ) can be calculated using 



I n (r,t) ex \ef\G 



; (UI 2 n 2 K 



t h c cot 0). (37) 



We have derived in this section general solutions de- 
scribing the spatiotemporal response of crystals in Bragg 
diffraction. In each particular case it is important 
to know the appropriate plane-wave response functions 
G 0H (^ H ). They can be calculated numerically in the 
general case, and examples are discussed in the follow- 
ing Sec. |III| for the reflection (Bragg) geometry, and in 
Sec. IV for the transmission (Laue) geometry. In some 



cases G ou (£ H ) can be calculated analytically, in particu- 



lar, for non-absorbing crystals. In Sees. |III| and |IV| we will 
derive analytical expressions for the response functions of 
non-absorbing crystals in the general case of asymmetric 
diffraction rj =/= 0, in Bragg and Laue scattering geome- 
tries, respectively, and perform analysis of the spatiotem- 
poral crystal response using these analytical solutions. 



III. RESPONSE IN BRAGG-CASE GEOMETRY 



A. Diffraction and Forward Diffraction Amplitudes 



We begin this section by summarizing the well-known 
results of the dynamical theory of x-ray Bragg diffraction 
for both the forward diffraction R 00 and diffraction R 0H 
amplitudes measured at the rear (z = d) and the front 
(z = 0) surfaces of a crystal, respectively: 



R nn =e 



XH i? 2 - Rl 



R 2 -i? lC iK -* a )d' 
1 _ e i(x 1 - x 2 )d 

R 2 -R^i-^W 



(38) 



where 



y= K ^[ba + Xo (l-b)}, 



(39) 



j\ = d/ A H . G = J\b\x H xJX B 



and 

A, 



Vi \ih\ 

sin 9 



A':, 



Aw 



sin 8 



K a \p\ vra^ 



(40) 



Here Xn ffl- = 0, H, H) are Fourier coefficients of the 
periodic-in-space crystal electric susceptibility xi r )- in 
general, x n are verv small complex parameters. The 
imaginary part 9{x M } is related to the cross-section 
of photo-absorption, while the real part is pri- 

marily related to the atomic Thomson scattering am- 
plitude. In many interesting cases, e.g., for Si crystals, 
3?{x H } 3> 3{x M }- In certain cases, e.g., for diamond 
or Be crystals, one can even neglect photo-absorption to 
a certain extent, and assume Xu t° be purely real pa- 
rameters. We make this approximation in the analytic 
calculations of the response functions below, although as 
shown by Kato [23] the resulting expressions can often be 
applied to absorbing perfect crystals by letting the x H be 
complex if the appropriate branches of square roots, etc. 
are taken. Typically, 3?{x w } ~ 10~ 4 — I0~ 7 for Si and 
diamond crystals for 5-20 keV x-rays. 

The index v — 1,2, identifies two possible solutions for 
the correction x v of the in-crystal wave vector k = K + 
kz with respect to the vacuum wavevector K . 

The diffraction R 0H and forward diffraction R QQ ampli- 



tudes in ( 38 1 are essentially functions of one main param- 



eter, the normalized deviation parameter y (39). It con 



tains all the information on the magnitude of the photon 



9 



frequency w , the direction of its momentum K a relative 
to diffraction vector H and to the internal surface nor- 
mal z, the asymmetry factor 6, and other information 
pertinent to scattering geometry. 



The parameter A H ( 40 1 in ( 39 1 is an extinction length 
[4"T] . In Eq. (40) we also define the extinction length 
in symmetric scattering geometry AjJ , for which the 
asymmetry angle r\ = 0, see Fig. [ija). In this case 
7 = — j H = sin 9, and A H — > A# . An important feature 
of the symmetric version A# is that it is invariant for 
a given Bragg reflection, being independent of the pho- 
ton frequency uj or incidence angle 9 to good accuracy 
for crystals with small photo-absorption; is deter- 
mined solely by the diffraction vector H . This can be 



seen from (40), by using the fact that \h ^ V-^o an< ^ 
sin 9 ~ 2H/ K . The extinction length A^ determines the 
characteristic interaction length in Bragg diffraction from 
the atomic planes with diffraction vector H . Along with 
the crystal thickness d, the extinction length is another 
characteristic measure of length in Bragg diffraction. 



With the help of the above Eqs. (39), and using expres- 



sion ( 13 ) for the deviation parameter a, the following 



relationship can be established between y and £1: 



fi = -sgn{6} 



y_ 



+ W„UJ, 



where 



T A = 



2A, 



c \l H \ 
7| S ) = 2AW/c 



= 7>> ^sintf, 



and 



(b-l) 
2b 



2shT 



(41) 

(42) 
(43) 

(44) 



Here, T A is the characteristic measure of time in Bragg 
diffraction associated with the extinction length. It can 
be directly compared to T d — AT A pi]), another charac- 



teristic measure of time in Bragg diffraction associated 
with the crystal thickness. We have introduced here also 
the Bragg reflection invariant time constant 7l s) ( 43 1 , 



which is associated with the Bragg reflection invariant 
extinction length in symmetric scattering geometry Ah . 
Typically, A' H ' « 1 — 50 /im (see Table [j] in Appendix [e|) , 
and therefore 7?' ~ 50 - 1000 fs. 



The parameter w H in (44) is a Bragg's law correction 
due to refraction at the vacuum-crystal interface |18j . 
Its magnitude w„ in symmetric diffraction (b = —1) is 
a Bragg reflection invariant, similar to the invariance of 
the extinction length A ( ^' . While the precise value of Wh 
depends on the Bragg reflection, in most cases it is very 
small (w ( h < 10~ 4 , see Table [j] in Appendix [e]). 

Far off the region of Bragg diffraction, where the devi- 
ation parameter |y| 1 or equivalently when a pi) is 



large, the diffraction signal is R 0H (oo) — (38), as ex- 
pected. On the contrary, 



R oo(°°) = c > c = ex P iX 



K d 
•2% 



(45) 



has a non-zero value that represents the diffraction- 
free transmission amplitude of the incident radiation 
with refraction and photo-absorption accounted for by 
C through the complex x ■ The actual forward diffrac- 
tion amplitude is therefore obtained by subtracting off 
the trivial y- independent amplitude C (45), 



D D 



c, 



(46) 



resulting in i? 00 (oo) = 0. Using the actual forward 



diffraction amplitude ( 46 1 , the forward Bragg diffraction 



response function ( 18 ) can be presented as the sum 



G 00 (U =CS(Q + G 00 (Q 



f!.-«e.ju, + n), («' 



so that the response is decomposed into the prompt 
diffraction- free transmission C(5(£ ), and the delayed ac- 
tual forward diffraction response function G 00 (£ ) . A sim- 
ilar approach has been used in [21] to deal with time de- 
pendence of forward resonant scattering from Mossbauer 
nuclei. 



B. Response Functions 



The response functions G n H (£„ ) and G 00 (£ ) are calcu- 
lated using Eqs. (18) and (47) with forward diffraction 
R 00 and diffraction R 0H amplitudes given by Eqs. (38) 



and (46) of the previous sections. We use Equation (41) 



to perform the integration over y instead of ft in the 
Fourier integrals. In reflection (Bragg) scattering geom- 
etry 7 H < 0, the asymmetry ratio b < 0, and there- 



fore the relationship (41) between f2 and y is actually 

The response functions G 0H (£ H ) and G 00 (£ ) can be cal- 
culated numerically in the general case. Figure [4] shows 
examples of such calculations. The left column of Fig. [4] 
shows examples of reflectivity spectra |-R 0H | 2 for crystals 
of different thicknesses d. Crystals of rather large thick- 
ness are considered d 3> A H . In the particular case of 
the H = (004) Bragg reflection in diamond crystal, the 
extinction length A ( h = 3.6 /im. The reflectivity spectra 
are nearly rectangular within \y\ < 1, and have an almost 
crystal thickness independent form and width AE H . The 
reflectivity is almost 100% . Such a high reflectivity is 
typical for diamond crystals due to low photo-absorption 
and a high Debye- Waller factor [24l [25] . The second col- 
umn in Fig. [4] shows results of calculations for the for- 
ward diffraction intensity spectra. They looks like inverse 
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FIG. 4: Spectral dependences (left three columns) of the Bragg diffraction (BD) intensity \R 0H (E)\ , the forward diffraction in- 
tensity \R 00 (E)\ 2 , the actual forward Bragg diffraction (FBD) intensity |7? 00 (_B)| 2 = \Roo{E) — 7?oo(oo)| 2 , and the corresponding 
temporal intensity dependences of the response functions \G 0U (i)| 2 (right column) in symmetric Bragg-case geometry. Numeric 
calculations use Eqs. ( |38[ )-(44 1, (13 1, and Eqs. (18l-(21l for the glancing angle of incidence 8 = 45° to the (004) reflecting 
atomic planes in diamond, with the asymmetry angle r\ = 0. The center of the Bragg reflection region (y — 0) corresponds to 
x-ray photon energy E c = 9.83 keV. The spectral and time dependences in diamond crystals of different thickness d — 0.05 mm, 
d — 0.1 mm, and d — 0.2 mm are shown in three different rows from bottom to top, respectively. 



diffraction spectra, because of the dominating contribu- 
tion of the trivial transmission in the diffraction-free re- 
gion \y\ 3> 1. The third column shows intensity spectra 
of the actual forward diffraction. The main contribution 
is outside the region of the total Bragg reflection. The 
spectral width AE Q ^> AE H and is crystal thickness de- 
pendent, varying linearly with d. 

The last (right) column shows the temporal dependences 
of the diffraction response function intensity \G 0H (£ H )\ 2 , 
and the actual forward diffraction response functions in- 
tensity |G 00 (£ )| 2 . In agreement with the behavior of 
the spectral dependences, \G 0H (^ H )\ 2 is approximately 
independent of the crystal thickness, while \G 0H {£, H )\ 2 
strongly depends on the crystal thickness d. The charac- 
teristic times of Bragg diffraction T A and forward Bragg 
diffraction T are indicated on the graphs. The temporal 
features at £ H = T d for the crystal with smallest thickness 
d = 0.05 mm represents diffracted and forward diffracted 
wavefields originating from the rear and front surface re- 
flections as schematically illustrated in Fig. |3j c.f. also 



numeric calculations in [T2J [18] . In thicker crystals these 
echo wavefields arrive at later times, which are outside 
the presented time range. Interestingly, the characteris- 
tic time of the Bragg diffraction response in the range 
£h > 7~ d changes from T A to T . This reflects the fact 
that the contribution to this signal comes from the same 
modes which contribute to forward diffraction, i.e., from 
the modes propagating through the whole crystal thick- 
ness, and not from those propagating only through the 
extinction length. 

More insight can be obtained from analytical solutions. 
The response functions can be calculated analytically in 
some specific cases, e.g., in the approximation of a non- 
absorbing (3 {x n } — 0) and thick crystal, for which d ^> 
A H , or equivalently A 3> 1. In this case, the diffraction 
and forward diffraction amplitudes can be approximated 
by 



Ran = G 



-y + 



for \y\ < 1 



-V + sgn {y} \/y 2 - I for \y\ > 1, 



(48) 
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c 



for \y\ < 1 



: A 



-II - 



sgn \y\ vT _ 1 , 

WV " ' for \y\ > 1. 

(49) 

The following observations were used to obtain Eqs. ( [48] )- 
p9| . In the region y ^ 1, both i?i < 1, and Ri/R 2 < 1, 
while for y ^ -1, both i? 2 < 1, and i? 2 /i?i < 1- 
Neglecting these small terms, Eqs. (38) transform to 
Eqs. (E^-pjl). 



Equation (48) represents a well known result of the 



dynamical theory, that Bragg diffraction from a non- 
absorbing, thick crystal i > 1 takes place with to- 
tal (100%) reflectivity \R 0H (y)\ 2 ^ 1 within the region 



\y\ < 1, or equivalently, using (41), within the photon 
energy range: 



(50) 



in agreement with results of numeric calculations shown 
in Fig.ga). Here AE H = hAfl. 

Using Eqs. (|46|-(49 ), the actual forward diffraction am- 
plitude can be presented in the |y| > 1 range as 



R„ 



cjcxp [i^( 



-y + sgn 



{y}V: 



y- 



l 



l 



Equation ( 51 ) is in agreement with the results of numeric 
calculations shown in the third column of Fig. [4j The 
forward diffraction spectral width AE ~ AE H (A/2ir) 
is a factor of (A/2tt) broader than the Bragg diffraction 
spectral width, and is crystal thickness dependent. 



In general, however, an asymmetric geometry changes 
the characteristic time constants. 



According to Eq. (52), the characteristic time constant 
in Bragg diffraction is T A (42), which is a function of 
the asymmetry factor b; in fact, it scales with i,y\b\. By 



appropriately choosing the asymmetry factor 6, the time 
response can be made faster or slower compared to the 
time response of Bragg diffraction in symmetric geome- 
try. Additionally, the uncertainty- type relationship ( 50 ) 



associates the characteristic time constant of diffraction 
with its spectral width. As a consequence of Eqs. (50) 



and ( 42 ) , a well know result of the dynamical theory can 



be reproduced: the Bragg reflection spectral width scales 
with l/vl&l, AE H = he / (A { h ,J\b\ sin 8). Note that both 
the typical energy and time scales are predominantly de- 
termined by a single parameter, namely, the extinction 
length AJ/ . This fact explains why Bragg diffraction is 
not instantaneous, as it builds by multiple, coherent scat- 
tering of x-rays within the extinction length. 

The characteristic time constant of forward Bragg diffrac- 
tion is significantly different, since different multiple scat- 
tering processes are involved. According to Eqs. (53)- 
( 54 ) , the characteristic time is given by T , which is a 



factor of A = d/A H smaller than the characteristic con- 
stant of Bragg diffraction T A ■ Interestingly, T is practi- 
cally the same as in symmetric scattering geometry. In 
the general case, it is basically defined by the Bragg re- 
flection invariant A« and the effective crystal thickness 
seen by incident x-rays d/~f . From Eq. (53), we also 



calculate that the first trailing maximum of the forward 
diffraction response function appears at t s — 267^,, as 
illustrated in Fig. [1J and its duration is At s — 16.57^,. 

f49|), |5l|), and relationship We note also that, if £ <C T d , the expression for the 



Using the diffraction amplitude R 0H (48), the forward 
diffraction amplitude R Q0 
Q = y/T A + w h lu (42), we obtain the plane- wave re- 
sponse functions Eq. ( 18 1 in the Bragg-case geometry (see 
Appendix [A] for mathematical details): 



forward diffraction response function given in Eqs. ( 53 ) 
( 54 1 can be simplified to 



(52) 



C Ji(VU% 

^ VUT 



(55) 



^00 (4 ) 



2T n 



(53) 



% = rjA 



c(d/l ) 



These solutions are valid if £ n < T d (H — 0, H), i.e., over 
the duration of time that is less than the total propa- 
gation time through the crystal T d - Thus, this solution 
does not include possible reflections from the rear and 
front crystal surfaces. 

In the limit of symmetric Bragg scattering, the response 
function envelopes are the same as those obtained in pQ. 



Analysis of the Spatiotemporal Response in 
Bragg-case Geometry 



By combining the analytical expressions for the plane- 
wave response functions (52)- (53 1 obtained in the previ- 



(54) ous section with the general solutions (36)-(37), we are 



now in a position to describe the spatiotemporal response 
of crystals in x-ray Bragg diffraction resulting from the 
excitation by an ultra-short and laterally confined x-ray 
pulse. 

To make the analysis more instructive, we show in Figs. [5] 
and [6] examples of 2D (r n , v n ) color plots of the spa- 
tiotemporal intensity profiles of forward Bragg diffraction 
(FBD) and Bragg diffraction (BD) from a 100 /mi thick 
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FIG. 5: Spatiotemporal intensity profiles of FBD and BD from a 100 /im thick diamond crystal from the (004) Bragg reflection 
(A ( ^' = 3.6 /jm), in the refl ecti on (Bragg) scattering geometry - Fig, [lj^a) - with asymmetry angles r\ = 20° (a), r\ = 0° (b), 
r\ = —20° (c). We plot Eq. (37 1 using the plane-wave response functions ( |52[ )-( 53 1 and a Gaussian lateral spatial profile of the 
incident x-ray beam with a v — 1000 fim (i.e., a practically unbounded incident wavefront), were used in the calculations. Other 
parameters: 8 = 45°, E = 9.8 keV, which are the same as those used for the calculations of the response functions shown in 
Fig. [4] The intensity front tangent in Bragg diffraction (BD) is dv H /dr H = —c/D (22 1. White dashed lines are traces of the 
Bragg's law dispersion envelopes H 2 (v n — r n ccot8) |36l-|37l with a tangent Av rL /Ar u — ccot#, c.f. Pig. |g] 



diamond crystal in the (004) Bragg reflection (A ( „ = 
3.6 fj,m) with the asymmetry angle i] = 20° (a), r\ = 0° 
(b), or 7} = -20° (c). 

Fig. [5] shows examples of calculations that apply the 
above mentioned equations to an incident wavefront that 
is, for all practical purposes, laterally unbounded (we 
assume that the incident spatial profile has a Gaussian 
distribution with a v = 1000 fim). In the symmetric case 
- Fig.[5jb) - the spatiotemporal profiles of both FBD and 
BD are homogeneous in the lateral spatial shift v H , i.e., 
they show no variation along the plane perpendicular to 
the appropriate optical axis u H - see Fig. Qa). The FBD 
response remains independent of v H for non-zero values of 



the asymmetry angle 77 7^ 0, as the plots demonstrate in 
the left columns of Figs.[5][a) and (c). In contrast, the BD 
profiles acquire modulations along v H if rj 7^ 0. They also 
produce the impression that the wavefronts of the BD 
wavefields are inclined. There are two phase factors in 
the expression for the wavefield £ H (r,t) in Eq. ( 36 1 . The 
first is exp(— icjr K ) defining the wavefront perpendicular 
to the optical axis u H , and another one exp(— iw H ui£ H ) 
resulting from the plane- wave response function G 0H (£ H ) 
( 52 1 . Since the second contribution is due to a small 



refractive correction, the wavefront is practically not in- 
clined. The pronounced effect seen in Figs.[5|a) and (c) is 
actually the inclined amplitude (intensity) front due to 
the amplitude modulation perpendicular to u H result- 
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FIG. 6: Spatiotemporal intensity profiles of FBD and BD for similar parameters as that of Fig. [5] Here, however, the incident 
x-rays have a bounded lateral spatial profile Tl(v ) that is assumed to be Gaussian with a v = 10 /im. Due to this, the 
spatiotemporal intensity profiles shown in Fig. [5] are limited now by a tight Bragg's law dispersion envelope II 2 [v H —t^c cot 8) 
(37 1 with a tangent dv n /dT n = ccot#, shown as white dashed line. 



ing from angular dispersion, discussed in Sec. |IIB| For- 
mally, the inclination and modulation reveal themselves 
through the argument £ H of G 0H (£ H ) which depends both 
on time t, and, if 77 7^ 0, also on the space variable v H - 
Eq. (23). The magnitude of the inclination to u H is D/c, 
it scales with the normalized angular dispersion rate D. 
The inclination of the intensity front changes sign with 
the sign of r\. The tilting of the intensity profiles due to 
Bragg diffraction was previously noted by Bushuev |20j . 

In all cases, varying the magnitude and sign of the asym- 
metry angle r\ changes the time constants T and T A , 
resulting in either dilation - Fig. [5]ja) - or contraction 
- Fig. [5^c) - of the oscillating intensity structures as- 
sociated with the spatiotemporal response; this dilation 
or contraction as compared to the symmetric case oc- 
curs both in the w„ and the t\, directions, as shown in 



Fig.§b). 

In the next step, we narrow considerably the lateral spa- 
tial profile of the incident x-ray beam. Figure [6] shows ex- 
amples of calculations for incident x-rays having a Gaus- 
sian lateral spatial profile with <j v = 10 fim, with all other 
parameters being identical to those in Fig. [5] The wave- 
field in the direction perpendicular to the optical axis 
u H is bounded by the Bragg's law dispersion envelope 
U 2 (v n ~ T n ccot0) - Eqs. @-([37|. White dashed lines 
in Figs. [3] and [6] are traces of the envelope. The tangent 
dv n /dr n = c cot 9 is independent on whether the geome- 
try is symmetric 77 = - Fig. [6] (b) , or asymmetric rj =/= 
- Fig. [6] (a), (c). As has been mentioned in Sec. II C this is 
a result of Bragg's law dispersion due to angular spread 
in the incident beam caused by the bounded wavefront. 
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The two effects of lateral amplitude modulation of the 
wavefield, due both to angular dispersion and to Bragg's 
law of dispersion, can be clearly distinguished by com- 
paring the spatiotemporal profiles in Fig.[5]and in Fig. [6] 



IV. RESPONSE IN LAUE-CASE GEOMETRY 



Here we use characteristic time constants defined pre- 
viously^ T d in Eq. pTJ, T A in Eq. ggh, and % in 



fied if f H <T d : 



C 00 (£ ) 



Eq. (54). Both equations (59) and (60) can be simpli- 



2 % ^T 



(61) 



A. Diffraction, Forward Diffraction Amplitudes 
and Response Functions 



-iw H w£ H 



(62) 



The wavefield amplitudes in transmission (Laue) geome- 
try are given by the following expressions [Still]: 



Rn, 



1 



R-2 — R\ 
R1R2 



R 2 e lx i d 



R x e l 



Ro — R 



e i>r i rf 



(56) 



The notation is the same as in ( 39 ) , but in contrast to the 



reflection (Bragg) geometry, the asymmetry factor b ( 39 ) 
is positive in transmission geometry. Using Eqs. (39)- 



( 44 ) , the wavefield amplitudes ( 56 1 can be presented as 
R 00 = C e-^f/ 2 W(y), 



W(y) = cos 



A 



sm 



m ■ 



R 0H = iC Q e 



-LAy/2 



V(y), 



(57) 



sm 



V(y) = 



{iVi 



(58) 



Using the forward diffraction amplitude R 00 given by 
(57), the diffraction amplitude R 0H given by (58), the 



relationship SI = —y/T A + w h lu (42), and assuming zero 



photoabsorption, we compute the plane-wave response 
functions Eq. ( 18 ) for the Laue-case which are given by 
(see Appendix Bjfor mathematical details): 

Goo (4) =G 00 (Q +C6&), 

J, 



[o<£o <r d ], 



. C Q 

G ff (Cff ) = _ 1 ^0 



1 - 



M 1 ® 



(59) 




1 - 



T, 



[o < e„ < rj. 



(60) 



Remarkably, for small < 7^, the forward diffraction 
plane- wave response function G* nn (g n ) in the Laue-case 



geometry (62), and its counterpart (55 



geometry, are equivalent, however, wit 
a consequence of sgn {b} in Eq. (41 ). 



in the Bragg-case 
1 inverted signs, as 



The time constant 7~ ( 54 ) is essentially the same in both 



transmission (Laue) and reflection (Bragg) geometries. It 
equals the time constant 7l B) in symmetric Bragg diffrac- 
tion scaled by a ratio Ah /(d/~f ) of the symmetric extinc- 
tion length Ah to the effective crystal thickness (d/%), 
i.e., the crystal thickness seen by incident x-rays. The 
time constant of forward Bragg diffraction T ( 54 ) is thus 
general for all symmetric or asymmetric, transmission or 
reflection scattering geometries. The primary parameter 
controlling the forward Bragg diffraction response is the 
effective crystal thickness. 



B. Analysis of the Spatiotemporal Response 



Expressions (59)-(60) for the plane-wave response func- 



tions are very similar to the analogous expressions ob- 
tained by Graeff and Malgrange in [T5] [TB] , with the ex- 
ception that G 00 (£ ) in (59) contains also the prompt 
(5-function contribution. It originates from the spec- 
tral components in the incident pulse with frequencies 
far from the Bragg diffraction region that propagate 
diffraction-free through the crystal. In Appendix [B] wc 
also provide some more details on the comparison with 
the results of [T31[T5]. 

Figure [7] shows results of numeric calculations of the 
plane- wave response function intensities \G on (£_ H )\ 2 and 
-iw H uj£ related to them spectral dependences |i? 0K (-E)| 2 for the 
Laue case. These calculations are provided to facili- 
tate "visualization" of the analytical solutions given by 
Eqs. (56)-([60|). The dependences shown in Fig. 7] are 



counterparts of the analogous Bragg-case dependences 
shown in Fig. [4j They are calculated under the same con- 
ditions, with the single difference being that the asym- 
metry angle is now r; = 7r/2 ("symmetric" Laue case), 
instead of 77 = in Fig. [4] (symmetric Bragg case). 

In the Laue-case, the spectral range where Bragg diffrac- 
tion takes place scales with AE H . While this is similar 
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E-E c [eY] E-E c [eV] E — E c [eV] 



FIG. 7: 



Spectral dependences (left three columns) of the Bragg diffraction (BD) intensity \R 0H (E)\^ 



the forward Bragg 

diffraction intensity |i? 00 (J5)| 2 , the actual forward Bragg diffraction (FBD) intensity |_R 00 (_E)| 2 = ji?oo(£') — J?oo(oo)| 2 , and 
the corresponding temporal intensity dependences of the response functions \G 0U (£ n )\ 2 (right column) in Laue-case geometry. 
Numeric calculations use Eqs. (561, ( 39 (- ( 44 1 and ( 18 1 for the glancing angle of incidence 6 — 45° to the (004) reflecting atomic 
planes in diamond, with the asymmetry angle r\ = 90°. The center of the Bragg reflection region (y = 0) corresponds to the 
x-ray photon energy E c = 9.83 keV. The spectral and time dependences in diamond crystals of different thickness d = 67 /im, 
d — 80 (jm, and d = 104 /im are shown in three different rows from bottom to top, respectively. 



to the Bragg-case geometry, for Laue there is no region 
of total reflection. The dominant feature of the spec- 
tral intensity dependence in Fig. [7] is the intensity os- 
cillations, which is associated with the well-known Pen- 
dellosung effect of Ewald |4]. Pendellosung is basically 
related to oscillations of the energy flow between the for- 
ward diffracted and diffracted beams in the Laue-case 
geometry. The period of oscillations scales with the ex- 
tinction length A H , which is 3.6 /xm in this particular 
case. For the calculations presented in the middle row 
the crystal thickness d is reduced by 2ttA h = 22.6 /j,m, 
compared to the d value used in the calculations pre- 
sented in the top row. The crystal thickness d is fur- 
ther decreased by ttA h = 11.3 /xm for the calculations 
presented in the bottom row. Varying the crystal thick- 
ness leads to periodic in d and in E variations of the 
intensities that are complementary for the diffracted and 
forward diffracted signals. Interestingly, in the spectral 
range \E — E c \ > AE H , the actual forward diffraction 
intensity \R 00 (E)\ 2 = \R 00 (E) — -Roo(°°)| 2 has a struc- 
ture very similar to that in the Bragg-case geometry, c.f., 



Fig. 7J This is in agreement with the fact that G 00 (£ ) for 
small £ values are identical (modulo the inverted sign) in 
Bragg-case and Laue-case geometries, c.f. Eq. (55) and 
Eq. (611. 



Comparison of Eq. ([60]), and Eq. (|52j), as well as the 
results of numeric calculations of \G 0H (£ H )| 2 in Figs. [7J 
and |4j show that the characteristic time of diffraction 
in Laue-case geometry is T (54), i.e., different from the 
characteristic time of diffraction T A ( 43 ) in Bragg-case ge- 



ometry. This evidences that two different characteristic 
length scales are involved for these two different diffrac- 
tion cases. 



Eqs. ( 59 )-( 60 ) and Fig. [7] demonstrate a signature feature 
of the Laue-case plane- wave response functions G on (£ M ) . 
Unlike the Bragg-case analogs, G (m (£, n ) vanish outside 
the range > £ n > T d . This effect has been reported 
and discussed by Shastri et al. [12l [13] using numeric 
calculations, and by Graeff and Malgrange using analyt- 
ical solutions in [T5J [TB] ■ 




400 600 

= t — u H /c [is] 

FIG. 8: Spatiotemporal intensity profiles of FBD and BD from a 100 /im thick diamond crystal from the (004) Bragg reflection 
(A§ = 3.6 fim), in the transmission (Laue) scattering geometry - Fig.JfJb) - with asymmetry angles r\ = 110° (a), r\ = 90° (b), 
r\ = 70° (c). We plot Eq. (371 with the plane-wave response functions (|59|-( 60 1 and a Gaussian lateral spatial profile of the 



incident x-ray beam with a v = 1000 /im, (i.e., a practically unbounded incident wavefront) were used in the calculations. Other 
parameters: 6 = 45°, E = 9.8 keV, which are the same as those used for response functions shown in Fig. [7] The intensity 
front tangent in Bragg diffraction (BD) is dv n /dr n = —c/D (221. White dashed lines are traces of the Bragg's law dispersion 
envelopes H 2 (v u — T n ccot8) (361-1371 with a tangent dv n /dr n = ccotO, c.f. Fig. [9] 



This feature, however, deserves a more detailed discus- sion, as it is in fact valid only under certain conditions, 



17 




T Q =t- Ujc [fs] T H = t - ujc [fe] 

FIG. 9: Spatiotemporal intensity profiles of FBD and BD for similar parameters as that of Fig. [8] Here, however, the incident 
x-rays have a bounded lateral spatial profile n(w ) that is assumed to be Gaussian with a v = 10 /im. Due to this, the 
spatiotemporal intensity profiles shown in Fig. [8] are limited now by a tight Bragg's law dispersion envelope II 2 (v n —T n c cot 9) 
( |37| l with a tangent dv H /dr u = ccot(9, shown as white dashed line. 

but not in general. To illustrate this, we refer to the ex- poral response of forward Bragg diffraction (FBD) and 
ample 2D (r n , v n ) intensity color plots of the spatiotem- Bragg diffraction (BD) in Laue-case geometry shown in 
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Figs. M and |9l which are analogs the Bragg-case shown in 
Figs. [5] and [6[ 

Fig. [8] shows 2D plots for practically unbounded incident 
wave wavefront. Unlike the Bragg case, in the Laue-case 
geometry the intensity fronts of the diffracted wavefields 
are always strongly inclined. This is a consequence of 
the non-vanishing angular dispersion in Laue geometry, 
in agreement with Eqs. (22)-(23). 



Laue-case FBD is truly limited in time, both for obser- 
vations made at a single point or over an extended field 
of view: the duration is always T d (21 1, which can be 



changed by d, 9, and r\. The Laue-case BD is limited in 
time for an observer who measures the field at a single 
point, and its duration is the same as in FBD, namely, 
T d . However, for an observer that collects the reflected 
x-rays over some region in space, the duration of BD de- 
pends on the extent of the field of view. For an infinite 
field of view the duration is infinite. 

If the incident wavefront is now strongly bounded, as as- 
sumed for calculations of the 2D plots presented in Fig.[9j 
then BD is limited to the region of the Bragg's law dis- 
persion envelope. As a result, the duration of BD for an 
observer with an infinite field of view becomes limited 
to T d /b. We obtained this result also by ray tracing the 
wavefronts in Fig. [2] which represents the limiting case 
of an extremely bounded incident wavefront. 

From the above examples it is clear that there is no un- 
ambiguous answer to the question what is the duration 
of x-ray diffraction in Laue-case geometry. Depending 
on the conditions of the experiment, it can be either T d , 
or T d /b, or even arbitrarily long. The duration of BD in 
the Laue-case can be varied not only by decreasing the 
crystal thickness, as was suggested in [15], but also by 
varying T d or T d /b through the asymmetry angle 77 and 
asymmetry factor b, as follows from (21). 



The lateral spread v n of FBD and BD in Laue-case geom- 
etry appears to be limited to v^ ax ^ if the incident beam 
has a bounded wavefront. This is a well known result of 
the dynamical theory, supported by many experiments, 
reviewed in detail, e.g., in [TT]. Fig. [9] demonstrates how 
the limited-in-time crystal response correlates with the 
limited lateral spread. Using these graphs one can find 
that the maximal lateral spread is given by 



:J (rnax) = dcos2 6»/ 7ff , V ^ ax ) = dcOs26>/ 7o (63) 



The same values can be obtained using ray tracing of the 
wavefronts in Fig. [2] In agreement with ( |63[ ), we find 
from Fig. ^AA' = V ( max *> = dcos26/~y H , and BB' = 

v (vaax) _ dcOs26>/7 . 



V. APPLICATIONS 



A. Self-Seeding of XFELs 



Understanding spatiotemporal dependencies in Bragg 
diffraction of x-rays has immediate practical implica- 
tions, in particular for self-seeding of x-ray free-electron 
lasers (XFELs). The self-seeding scheme uses an up- 
stream XFEL to generate an intense x-ray pulse via self- 
amplified spontaneous emission (SASE). The relatively 
broad-bandwidth SASE pulse is then put through an x- 
ray monochromator to generate a monochromatic seed 
for the downstream XFEL undulators, which in turn am- 
plifies the narrow bandwidth seed to produce fully co- 
herent x-rays [261 12 7j . However, traditional two- or four- 
bounce monochromators induce a large delay (> 10 ps) 
of the x-rays, which in turn requires an impractically long 
(~ 40 m) electron beam transfer line. 

A very clever, readily realizable idea of a "wake" 
monochromator that produces a monochromatic x-ray 
seed at an optimal w 20 fs delay has been proposed by 
Geloni et al. [28, 29], and recently realized at the LCLS 
XFEL by an international team lead by Emma [30]. 

In the original proposal j28[ [29] the authors applied the 
equations of the dynamical theory of x-ray diffraction 
in crystals to calculate numerically the time dependence 
and strength of the monochromatic seed propagating in 
forward direction. The action of the monochromator 
crystal in the Bragg-transmission geometry was inter- 
preted in terms of a Bragg diffraction (BD) band-stop 
filter. The underlying physics is actually related to for- 
ward Bragg diffraction (FBD). We have discussed in de- 
tail its properties in the symmetric Bragg-case geometry 
relevant for self-seeding in [T]. We showed that, first, the 
characteristic time for FBD is T a (54 1, substantially dif- 



ferent (shorter) than the characteristic BD time T A ( 43 1 , 
and therefore the crystal in FBD generates a seed with a 
broader spectrum than a BD band-stop filter would do. 
Second, it was shown that the intensity of the monochro- 
matic seed is oc 1/T 2 , c.f. Eq. (55), which can there- 



fore be enhanced by varying parameters composing T Q - 
Similarly, its time delay t s = 26T (see Fig. [4j , and its 
duration At a = 16.57^ can be tailored by changing T , 
which can be done practically by adjusting the extinction 
length A ( „ (for example, by choosing another reflection or 
asymmetry parameter), or by changing the crystal thick- 
ness. A limitation of this scheme has been also identified 
in pp. It is due to the lateral shift of the FBD signal. 
This is a very generic effect, caused by the Bragg's law 
dispersion, as discussed in Sec. II C| of the present paper. 



The theory developed in the present paper allows us to 
diversify the variety of possible forward diffraction self- 
seeding monochromator schemes. First of all, forward 
Bragg diffraction in Laue-case geometry is a competitive 
approach. The possibility of applying FBD in Laue-case 
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geometry for self-seeding becomes immediately apparent 
from the derived equivalence of the forward diffraction 
plane- wave response functions \G 00 (£ )| 2 in the Laue-case 
geometry (62 1, and of its counterpart (551 in the Bragg- 



Position Sensitive 
Detector 



case geometry. The equivalence holds for small £ H -C T d , 
which is the range most appropriate for self-seeding of 
femtosecond long XFEL pulses. As has been established 
in this paper, the time constant of forward Bragg diffrac- 
tion T (54 1 is common for all symmetric or asymmet- 



ric, transmission or reflection scattering geometries, and 
is the only parameter which defines the strength, delay, 
and duration of FBD and therefore of the monochromatic 
seed. These properties advance FBD both in Bragg and 
Laue-case geometries, including asymmetric ones, to a 
universal approach for the generation of monochromatic, 
delayed seeds for self-seeded XFELs. The physics is con- 
trolled by the parameters which compose T (54): the 
magnitude of the effective crystal thickness dpy a , and 
the extinction length A ( ^' in the symmetric Bragg reflec- 
tion. Table [I] in Appendix [E] provides some useful data 
for Bragg reflections in diamond, which can be used to se- 
lect the Bragg reflection most appropriate for the desired 
application. Similar data for silicon and AI2O3 crystals 
can be found in |18j . 

There is no universal answer to the question: which ge- 
ometry is better, Bragg or Laue? We investigate certain 
aspects of this question below, where for simplicity we 
have restricted our analysis to the symmetric diffraction 
geometries, denned by r\ = in the Bragg-case - Fig.JTJa), 
and by 77 = 7r/2 in the Laue-case - Fig. [IJb). 

If the experimenter highly values operating the self- 
seeding monochromator over as large a spectral tuning 
range as is possible, than the Laue-case geometry may 
be a better choice. The strongest variation of the photon 
energy E with the glancing angle of incidence to the re- 
flecting atomic planes 8 takes place for small 8 < 7r/6, i.e., 
in the linear range of Bragg's law E sin = E H m. We 
write here Bragg's law in terms of photon energy E and 
Bragg energy E H = Hch/2, the smallest photon energy 
for which Bragg's law can be fulfilled (at 6 = 7r/2). In 
the symmetric Laue-case geometry the effective thickness 
d/-f = dj cos 8 does not vary much if is small, unlike 
the Bragg-case in which d/-f a = d/ sn\0. Therefore, a 
large variation in E is accompanied in the Laue-case ge- 
ometry with a small variation in d/j Q and therefore in 
7~ , resulting in a rather stable seed power and time de- 
lay of the seed over a large range of photon energies. In 
the Bragg-case geometry this is not the case. In addi- 
tion, the Laue-case geometry at small 8 allows for using 
thicker crystals for the same T a as compared to the sim- 
ilar situation in the Bragg-case. This may represent a 
technical advantage since the fabrication of thin crystals 
is typically more challenging. 

Using small 8 angles, however, also has its disadvan- 
tages. The lateral spatial shift, given by the Bragg's 
law dispersion envelope LI 2 (v n — r n c cot 8) - Eq. ( 37 1 , 




FIG. 10: Schematic of ultra-fast time measurements by map- 
ping time delay on space in asymmetric x-ray Bragg diffrac- 
tion from a crystal. The intensity front is rotated by an angle 
/3 upon asymmetric Bragg reflection of x-rays. Here tan /3-D 
is determined by normalized dispersion rate D (22 1. 



is v H = T H ccot8, i.e., proportional to cot 8, and is max- 
imal in the range of small 8. This may not be signifi- 
cant for very short x-ray pulses that can use short delay 
times t s . However, if one wants to seed long XFEL pulses 
~ 50 — 100 fs, then the Bragg-case scattering geometry 
close to backscattering 8 — > ir/2 would be a more advan- 
tageous option, albeit at a decrease in the spectral tuning 
range. 



B. Ultra- fast Time Measurements by Mapping 
Time on Space 



Angular dispersion in asymmetric Bragg diffraction re- 
sults in an inclined intensity front of the diffracted wave- 
fields. This effect is illustrated in Figs.[5]and|8j in Bragg- 
case and Laue-case geometries, respectively. Inclination 
of the intensity front in asymmetric x-ray diffraction ge- 
ometry was, to our knowledge, first explicitly derived us- 
ing ray tracing in [31], where it was proposed to be used 
for x-ray pulse compression. 

Here, we suggest using the effect of intensity front in- 
clination for time measurements of ultra-fast processes. 
The schematic drawing in Fig. |10| explains the idea. The 
inclined intensity front allows mapping time onto space, 
as different parts of the inclined front traverse and thus 
probe the sample at different times. 

We assume that the process under study is triggered ho- 
mogeneously over its extent by, for example, an external 
laser. If we then measure the sample with an inclined 
x-ray intensity front from an asymmetric crystal, than 
different transverse positions will be probed at different 
times, so that time dynamics can be extracted with a 
spatially resolved x-ray detector. We assume that the x- 
ray wavefront is sufficiently broad so that we may ignore 
any small additional spatial shifts that take place due to 
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Bragg's law dispersion. 

To estimate the achievable time resolution, we first ne- 
glect the finite duration of the incident x-ray pulse, the 
finite duration of Bragg diffraction, and the sample thick- 
ness. In this case, the resolution of the time measure- 
ments At = (Ax/c)/ tan/3 = Ax/(cD) is determined by 
the spatial resolution of the detector Ax and normalized 
dispersion rate D (22 1. With Ax ~ 1—10 fim, and D ~ 1, 
an estimate for the time resolution is At ~ 3 — 30 fs. 
The duration of the incident x-ray pulse, the duration of 
Bragg diffraction, and the sample thickness will increase 
this number. The duration of the Bragg diffraction is 
T A in Bragg-case or T d in Laue-case geometry. By an 
appropriate choice of A„ and asymmetry factor &, the 
characteristic time Bragg diffraction (43) can be made, 
however, as small as T A — 1 fs, i.e. smaller than the At. 
Tilting the sample to the x-rays propagation direction 
(without tilting the intensity front) may result in a sim- 
ilar effect. However, tilting the intensity front in many 
cases may be advantageous, as it is decoupled from the 
propagation direction in the sample, which may be an 
important parameter of experiments, of diffraction ex- 
periments in particular. 



VI. CONCLUSIONS 



of Bragg diffraction in asymmetric geometries. Bragg's 
law dispersion produces a lateral spatial shift v a of pho- 
tons emerging from the crystal with respect to the inci- 
dent x-ray pulse that increases linearly with time delay 
r . A simple general relationship cr — v tan 8 holds in 
all diffraction geometries. This effect can be interpreted 
also in terms of the energy flow of the wavefields in the 
crystal. 

The spatiotemporal plane-wave response functions in 
Bragg diffraction can be expressed in terms of three char- 
acteristic space and related to them three time parame- 
ters: (i) crystal thickness d and T d (21 1, (ii) extinction 
length A H (40) and T A (43), (iii) rescattering length A^d 
and T (54). The glancing angle of incidence 9 and the 



asymmetry angle r\ also enter the three time parameters, 
and, therefore, are factors that change the spatiotempo- 
ral response scale. 

We address some practical applications of the developed 
theory. We show that forward Bragg diffraction (FBD) of 
x-rays in Laue-case geometry can be used for self-seeding 
of hard x-ray free electron lasers, along with FBD in the 
Bragg-case geometry. Laue-case FBD is advantageous if 
a large spectral tuning range is required. We discuss also 
a possibility of using asymmetric diffraction for ultra-fast 
time measurements with femtosecond resolution. 



The spatiotemporal response of crystals in x-ray Bragg 
diffraction resulting from excitation by an ultra-short, 
laterally confined x-ray pulse was studied theoretically. 
The theory developed in the paper presents an extension 
of the analysis in symmetric reflection geometry [1] to 
a generic case, which includes Bragg diffraction both in 
reflection (Bragg) and transmission (Laue) asymmetric 
scattering geometries. 

The spatiotemporal response is presented as a product 
of a crystal-intrinsic plane wave spatiotemporal response 
function and an envelope function defined by the crystal- 
independent transverse profile of the incident beam and 
the scattering geometry. The diffracted wavefields ex- 
hibit amplitude (or intensity) modulation perpendicular 
to the propagation direction due to two effects: angular 
dispersion and dispersion due to Bragg's law. Angular 
dispersion results in the inclination of the intensity front 
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Appendix A: Response Function in Reflection 
(Bragg) Geometry 



To compute the temporal response of the forward 
diffracted wave, we must evaluate 



G 00 (£ ) = ~^ e " 

' A 



2?r 



R 00 (y) 



(Al) 



where the R 00 (y) is given by subtracting C from (49). Since causality requires G 00 (£ < 0) = 0, we have found that 
the most convenient way to treat this particular problem is as an inverse Laplace transform. In the table of inverse 
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Laplace transforms given by Erdelyi, Magnus, Oberhettinger, and Tricomi [35], we find that 

ab 



dt e- pt f(t) = 1 - e -k(\/p 2 + Q2 -p) 



fit) 



o 



y/t{t + 2b) 



ay/tit + 2b) 



(A2) 



Now, we make the replacements: t — £ , p = —iy/T A , a = l/T A , and b — AT A /2; then, (A2| is proportional to the 
approximate Bragg transmission. Thus, we have that 



A J i [VU-AT A +Q/T A 
2" VUAT A +Q~ 



In terms of the forward Bragg diffraction amplitude, (A3 1 implies that 



2tt 



y ^oo(y) 



_c_ 
27: 



J, 



^ fi , - 



In 



with 7J, = T A /A = 2"f [A l S] 2 /(cd). On the other hand, the reflected wave is given by 



G h(£h) = t^c 

I A 



2tt 



'* )y R 0H (y), 



and the integrals are simple enough for Mathematica to do; we find that 

-{€ H y/T A ) (y- \fy 



dy . 



: > _ ! ) = j ^^/^^ gn^ /T ) + j for Z 7 ^ ) cos (^ - sin ^n/T A ) 



2^„/r A 



and 



1 dy 

.i 2tt 



] ys™(£ H y/T A ) + i\/l - y 2 cos{£ H y/T A ) 
Adding these two, we obtain 



= /i(WT A ) _ ,(ur A ) cos(e H /r A ) - sm(^/r A ) 



2^/r A 



^/rj 2 



= f e- iw H-^ r ^ e -i (g g /TJ y ^ ( ) = i g e -to g ^H Jl e(g H /T A ) 
'a J-oo 2tt 7 a W'a 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 



(A8) 



Appendix B: Response Function in Transmission and 
(Laue) Geometry 



To calculate response functions in Laue-case diffraction 
geometry we use Eqs. (18), (57), (58), and Eq. (41 ) in the 



form Q — —y/T A +w h uj, as b > 1 in Laue-case geometry. 
As a result we obtain 



G 00 (Z )=-^e-™*^°l 

' A 



x n = 



2ir 



W(y), Co = - 



A 

2 ' 



A 



W(y) = cos ( — \/y 2 + 1 ) + iy 



(Bl) 



G h\£h) 



-i — C^H^Hj 

q~ h j 

'a 

ge-iC^^y), C H 



^(y) =- 



sin 



Vy^TT) 



x/y^TT 



T A + 2' (B2) 



The Fourier integral 7^ in (B2) is a tabulated integral 
[3 3) and can be calculated analytically, as was previously 
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carried out in solving similar problems [TTJ [TH1 [53] : 

i H = i. /„ (V(-W-d) e | + Ch ) e | - C *) . 

(B3) 

Here 0() is the Heaviside unit step function whose value 
is zero for negative argument and one for positive argu- 



ment. 



The Fourier integral X in (Bl) can be calculated using 
the property 1 = dl H /d{A/2) - 8I H /dC, QH EES! M 
resulting in 



1 =1 



I = - 



A/2 + C 



j, U{a/2Y-q e - + Co e --Co 



'-4 



(B4) 



Here the <5-function appears as a result of differentiating the step functions. Finally with (B3| and (B4) and definition 
of £ and Q H in (B1)-(B2) we arrive at the following analytical expressions for the plane- wave response functions in 
Laue-case geometry 



G 00 (U = G„o(£o) +C5(Q, 



^oo (Co 



k 

T, 



r, 



(B5) 



1 - 



T, 



G 0H {i H ) 



.eg . 



n 




(B6) 



With one exception, these expressions agree with relevant 
expressions obtained by Malgrange and Graeff in [To] , 
where diffraction of short x-ray pulses with infinite wave- 
front in the asymmetric Laue-case was studied analyti- 
cally. Unlike the expression for forward diffraction pre- 
sented in [TB], Eq. (B5 ) contains the delta- function, which 



represents the prompt response in the forward diffraction 
due to spectral components far from the Bragg diffraction 
region that propagate essentially diffraction-free through 
the crystal. 

We note also, that a reference system (a/, z 1 ) was used in 
[TS] attached to the crystal rear surface. Unlike this, we 
are using in our treatment for each diffracted wavefield 
its own reference system (u n ,v n ). We are also using a 
different approach to calculate the vacuum wavevector of 
the diffracted wave ([7|-([To|. Due to these differences, the 
expressions for the spatiotemporal variables £ H ( 20 )- ( 23 1 , 
and similar variables in |16) - Eq. (28) - may appear at 
a first glance to be very different. However, our detailed 
comparison shows that they are actually identical. So 
mathematically our results and results of paper |16| for 
the delayed parts of the response functions are in agree- 
ment, except for the delta function in Eq. (B5|. 



Appendix C: Mapping Time on Lateral Space Shift 



The relationship v H = t h c cot 9 between the time delay 
t u and the spatial shift v H in Bragg diffraction, repre- 
senting the trace of the Bragg's law dispersion envelope 
n(i; H — t h c cot 9) in Eq. (36 1, can be derived alterna- 



tively by combining Ewald's concept of the crystal wave- 
field [4] with the concept of energy flow introduced by 
von Laue )34j . 

These concepts lead to the following picture of physi- 
cal processes involved in x-ray Bragg diffraction in crys- 

(b), 
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and 



tals. It is illustrated graphically in Figs, 
schematically presenting diffraction in the Bragg-case ge- 
ometry and in the Laue-case geometry, respectively. 

An incident monochromatic plane wave with wavevec- 
tor K a excites monochromatic waveficlds in the crystal 
given by Eq. , where each wavefield is associated with 
a tie point on one of the brunches of the dispersion sur- 
face, which we number below by v. The energy flow for 
each wavefield in a perfect crystal is given by the wave- 
field Poynting vector [23], which is directed along the 
normal to the dispersion surface taken at the tie point 
of the surface representing the field [35]. The Poynt- 
ing vector is parallel to the vector of group velocity V v . 
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FIG. 11: Schematic presentation of two-beam x-ray Bragg diffraction from a crystal (a) in the reflection (Bragg) scattering 
geometry and (b) in the transmission (Laue) scattering geometry. The direction of the energy flow by the crystal wavefield is 
indicated by the vector of the group velocity V u , at an angle <\> v to the direction of the incident wave. 



As a result, different monochromatic wavefields propa- 
gate, first, along different paths of different lengths, and 
second, with different group velocities V v . At the exit 
surface, the wavefield breaks up into independent plane 
wavefields, one with the wavevector K propagating in 
the direction u , and another with the wavevector K H 
propagating in the direction u H . One should note that 
the concept of the energy flow works with one substantial 
limitation: it is not applicable in the total reflection re- 
gion in Bragg-case geometry, as there is no propagating 
through the crystal wavefields in this case. It is appli- 
cable, however, outside the total reflection region. Due 
to this, the wave with the wavevector K„ is shown in 
Fig. [ill a) propagating along the line starting at point A. 

and (b) we calculate for the spatial 
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Using Figs 
shifts v and delays r. 



light in vacuum c, in agreement with Eq. ( C3 ) . Using the 
relationships OB = OA cos(f) vl OC = OA cos(2# — (j> v ), 
and Eq. ( C3 ) , we obtain 



t h c cot 9, U = (0,H). 



(C4) 



This relationship is valid in the general case of asymmet- 
ric diffraction, both for Bragg and Laue scattering ge- 
ometries. It maps temporal onto spatial scales in Bragg 
diffraction, in agreement with the Bragg's law dispersion 
envelope H(v u — r H c cot 9) in Eq. (36). Uncertainty re- 



lationships are always valid, and therefore Eq. (C4) are 



actually applied not for absolutely monochromatic waves, 
or waves localized in time and space, but rather for wave 
packets with certain spectral and momentum distribu- 
tions. 



v =AB = OAsm4> v , t = °^-—, (CI) 

v y c 

OA OC 

v H =AC = OAsm(26-<j>„), t h = — . (C2) 

V v c 

The magnitude of the group velocity can be given, as 
derived in Appendix dDj), by 



V v cos 9 

c cos(6* — <j> v ) 



(C3) 



We use in Eqs. (CI I the fact that propagation along OB 
(4> v = 0) or OC{£ 



Appendix D: Wavefield Group Velocity in the 
Crystal vs Propagation Angle 



The group velocity vector V„ is given by [TTJ [341 [35] : 

Vv = C ^L±%^ ; (Dl) 



i + Rl 



where R v is defined in Eq. ( 39 ) , and v numbers brunches 
of the dispersion surface. 

The absolute value of the group velocity V v can be calcu- 



29) takes place with the speed of lated by taking the magnitude of Eq. (Dl) and recalling 
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that u u H 



cos 29: 



y/l + 2Rl cos26» + i?4 



(D2) 



The angle <p v between the direction of the group velocity 
vector V v and optical axis u is determined from the 
vect or sc alar product cos0„ = u V v /V v , which using 
Eq. (Dl ) becomes 



COS( 



1 + Rl cos 29 



y/l + 2R 2 U cos 29 + Rt 



(D3) 



Combining (D2)-(D3| we obtain the following relation- 



ship between the magnitude of wavefield group velocity 
V„ and its direction <f> v 



Vu = _ 
c COS 



cos ( 



(D4) 



Equation (D4) gives a physically reasonable result. The 
direction of the wavefield propagation under the Bragg 
diffraction condition is at = 9, in which case the group 
velocity V v = c cos 9 is less than speed of light in vac- 
uum. Far from Bragg diffraction conditions <f> — or 
cj) = 29, resulting in a reasonable solution V v = c. 



Appendix E: Bragg Reflections in Diamond 



h k I 




A (B) 




AE H 
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[meV] 


1 
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3.01034 


1.09 


8.17 
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2 


2 





4.91561 


1.98 


3.04 


106. 


3 


1 


1 


5.76401 


3.74 


2.20 


56.0 


4 








6.95161 


3.63 


1.51 


60.6 


3 


3 


1 


7.57532 


5.89 


1.27 


35.8 


4 


2 


2 


8.51391 


5.03 


1.00 


44.5 


3 


3 


3 


9.03035 


7.83 


0.89 


27.3 


5 


1 


1 


9.03035 


7.83 


0.89 


27.3 


4 


4 





9.83108 


6.41 


0.75 


35.9 


5 


3 


1 


10.2815 


9.82 


0.69 


22.6 


6 


2 





10.9914 


7.87 


0.60 


29.2 


5 


3 


3 


11.3961 


11.9 


0.56 


19.1 


4 


4 


4 


12.0404 


9.44 


0.50 


25.0 


5 


1 


5 


12.4110 


14.2 


0.47 


16.4 


7 


1 


1 


12.4110 


14.2 


0.47 


16.4 



h k I 




A (B) 




AE H 




h k I 


E H 
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AE H 








[keV] H 


xlO -5 
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[keV] 


H 


xl0~ 5 
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6 


4 


2 


13.0051 
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0.43 
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7 


1 
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0.40 


14.3 
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TABLE I: Allowed Bragg reflections hkl in diamond crystals and their parameters relevant to the present studies: Bragg energy 
E H = hc/2d H , the extinction length A^ (40 1, the Bragg's law correction ( 44 1, the energy width AE H (at 8 = 7r/2). 
These parameters are calculated using interplanar distance d H = 3.56712(2) A in diamond crystals at T = 298 K [361 137] . a 
Debye Temperature of 2230 T [38] , and anomalous scattering factors from [391 40! . 
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